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Tantalizing cosmological and terrestrial evidence suggests the number of light neutrinos may be 
greater than 3, motivating a careful re-examination of cosmological bounds on extra light species. 
Big Bang Nucleosynthesis constrains the number of relativistic neutrino species present during nucle- 
osynthesis, N^p N , while measurements of the CMB angular power spectrum constrain the effective 
energy density in relativistic neutrinos at the time of matter-radiation equality, N^ AB . There are a 
number of scenarios where new sterile neutrino species may have different contributions to AiVjjP N 
and AA^ 13 , for masses that may be relevant to reconciling cosmological constraints with various 
terrestrial claims of neutrino oscillations. We consider a scenario with two sterile neutrinos and ex- 
plore whether partial thermalization of the sterile states can ease the tension between cosmological 
constraints on A r ™ N and terrestrial data. We then investigate the effect of a non-zero neutrino 
mass on their contribution to the radiation abundance, finding reductions in ATV^ 13 of more than 
5% for neutrinos with masses above 0.5 eV. While the effects we investigate here could play a role, 
we nevertheless find that two additional light sterile neutrinos species cannot fit all the data at the 
95% confidence level. 



I. INTRODUCTION 

Neutrinos are unique among the known elementary 
particles in that their properties have often been first, 
and in general more, constrained by astrophysical and 
cosmological limits than by direct laboratory measure- 
ments. Already in the 1970s cosmological probes gave, 
first, a constraint on the neutrino mass based on esti- 
mates of the density of non-relativistic matter in the 
universe pQ, and then a constraint on the number of 
light neutrino species based on estimates of primor- 
dial helium production during Big Bang Nucleosyn- 
thesis (BBN) [2-4 . Today, both of these probes have 
reached impressive sensitivity, and have started yield- 
ing some tantalizing suggestions of the possibility that 
extra neutrinos may be present in nature. 

The radiation abundance in neutrinos and beyond- 
standard-model relativistic species is usually expressed 
as the effective number of relativistic species, 



where pf = (77r 2 /120)(4/ll) 4 / 3 T 4 is the energy den- 
sity of one standard-model massless neutrino with a 
thermal distribution, p 7 is the energy density of pho- 
tons and yO re i is total energy density in relativistic par- 
ticles. In the standard model, by the time of BBN only 
the three known neutrino species contribute to p le \, re- 
sulting in N c g — 3.046 [5]- This is slightly larger than 
three due to reheating via e + e~ annihilation. 

Extra radiation beyond the standard model (the so 
called "dark" radiation) , would cause an excess (which 
we label AiV ff ) above the standard model value of N e g . 
Although adding an extra light fermion could con- 
tribute AN e s = 1, most generally N e g is non-integer 



and varies with time, and depends on the physics at 
play. Specifically, lepton asymmetries [51 |7j, particle 
decay [8] [9], partial thcrmalisation of new fermions 
[TO] ITT] , the effect of a new MeV-scale particle on the 
active neutrino temperature |12) , non-thermal produc- 
tion of dark matter 13 and heavy sterile neutrinos can 
all lead to contributions to N c g that are not integer 
and/or change with time. Therefore we can hope that 
probing AN c g precisely at different epochs - namely, 
during BBN and at the formation of the CMB - could 
discriminate between different models. 

Recent measurements of iV e ff have hinted at a value 
of N cS > 3.046 (AAcff > 0). Constraints on N cS 
can be derived from measurements of the primordial 
4 He mass fraction, Y,= i^f^. at BBN, T ~ 0.2 
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MeV. Izotov and Thuan [II] find Y p = 0.2565 ± 
0.0010(stat.)±0.0050(syst.), giving 7V c B ff BN = 3.68±g;|g 
or iV c BBN = 3.80±g;|g, each at 2cr, depending on the 
choice of the neutron lifetime, and assuming no lepton 
asymmetry. These are both more than la from the 
standard model value. Other recent estimates of Y p 
[T5] ITS] and various analyses of N e g at BBN, e.g. [TH - 
RU)] , give for the most part central values more than la 
above 3. 

CMB constraints on N c g are derived from measure- 
ments of the angular power spectrum at larger scales 
near the doppler peak, which can be used to constrain 
the redshift of matter-radiation equality (T ~ 1 eV). 
With independent measurements of the total matter 
abundance, this can also constrain the radiation abun- 
dance at matter-radiation equality. The South Pole 
Telescope suggests a somewhat high value, N^ S MB = 
3.71 ± 0.35 gU]. While evidence from WMAP 9 of 
an excess in N e g is now reduced somewhat from ear- 
lier data sets [2JJ, the actual value depends strongly 
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on which other observational data sets are combined 
with the WMAP data, with a new best-fit value of 
N^ 13 = 3.26 ± 0.35 when all data sets are combined. 

Interestingly, bounds from terrestrial searches for 
new physics on the masses and couplings of new parti- 
cles invariably result in constraints on their contribu- 
tion to the cosmological radiation, providing indirect 
constraints on AN e g. Of particular interest are the 
recent hints of a fourth, sterile neutrino species from 
reactor neutrino experiments [221 [23] > calibration data 
from Gallium-based solar neutrino detectors [24H26] . 
and the Short Baseline (SBL) neutrino beam experi- 
ments LSND [57] and MiniBooNE [2"5H3"T] which search 
for — > D e and — \ v e oscillations. All these gener- 
ally support the existence of least one sterile neutrino 
with mass ~ 0.1 — 1 eV. This neutrino would be pop- 
ulated in the early universe via an interplay of oscilla- 
tions and scattering, thus increasing iV e ff- 

Although the possibility of extra radiation due to 
sterile neutrinos seems to be substantiated at the gen- 
eral level, detailed analyses of the data reveal tensions 
between datasets and leave open the question of what 
scenario is most favored overall. MiniBooNE observes 
a difference between the muon neutrino and antineu- 
trino disappearance rates, hinting at CP violating ef- 
fects. The latest measurements by MiniBooNE show 
less tension between their neutrino/antineutrino re- 
sults, although the 3 + 2 scenario still provides a better 
fit to the data [3T] . The simplest explanation for this is 
the existence of two sterile neutrinos families, although 
other data does not easily accommodate that possibil- 
ity and the improvement in the global fit to the data 
may simply be due to the additional free parameters 
in a 3 + 2 model [3"2Tl36| . Combined fits of SBL and 
cosmological data have found further tension when the 
sterile neutrinos are fully thermalized, with the level of 
tension depending on exactly which datasets are con- 
sidered [3"?| . 

Whilst the cosmological data appear to rule out 
N e s = 5, multiple sterile neutrinos can still be accom- 
modated if one or more of them are not fully thermal- 
ized at the time of BBN. The degree of thermaliza- 
tion depends on the sterile neutrino mass and mix- 
ing parameters, as constrained by SBL data. Neu- 
trino density evolution and partial thermalization in 
a 3 + 2 scheme has been studied by Melchiorri et al. 
[llj . finding tension both between the various terres- 
trial data sets, and between terrestrial and cosmolog- 
ical data. Since then, there have been substantial im- 
provements in cosmological measurements and experi- 
mental results. 

Considering that there is evidence from multiple 
sources that perhaps additional light neutrinos exist, 
and scenarios with two sterile neutrinos have been pro- 
posed as a way to explain the MiniBooNE results, it is 
important to re-examine cosmological constraints with 



a more careful eye. With this goal in mind we have ex- 
plored both partial thermalization of the sterile species 
at BBN and the effects that small neutrino masses will 
have on the relativistic neutrino fraction at the time 
of matter-radiation equality. As we will show, N e g at 
BBN and CMB can be quite different for light neutri- 
nos, but even incorporating this fact, and pushing all 
constraints to their 2a level, cosmology can still not 
accommodate models with two such neutrinos. 



II. PARTIAL THERMALIZATION IN A 3 + 2 
SCENARIO AND BBN CONSTRAINTS 

We consider a scenario with two sterile neutrinos, 
and study the sensitivity of their thermalization effi- 
ciency to their masses and mixings. Specifically, we 
denote the two sterile neutrino flavors as v s and i/ r , 
and the corresponding mass eigenstates as 7714, 7715, 
such that the mixing matrix U has entries U S 4 ~ 
U r s ~ 1, and the hierarchy m 5 > m 4 3> rrij (J 
1.2. ,' ! : holds. For simplicity, we also assume £/ r 4 = 
U T 5 — 0, so that our results for A^ BN only depend 
on J7 e 4, t/ M 4, U e s, Ufjs. We have verified that the com- 
plex phase 77 = Arg([/* 4 U e5 U* 5 ) does not affect 
the degree of thermalization, therefore we simplify the 
notation by considering U to be real. 

In the density matrix formalism the differential equa- 
tions governing evolution of the neutrino density are 



p = Hp- P Hl =i[H m + V oB , p]-{^,(p-p eq )}, (2) 

where p is the 5x5 neutrino density matrix in the flavor 
basis with diagonal entries corresponding to physical 
densities, H is the full Hamiltonian, H m = U H W 
is a rotation of the free neutrino Hamiltonian in the 
mass basis H a = diag(£'i, E 2 , E 3 , E4, E 5 ), and p eq 
is the density matrix at thermal equilibrium, p eq = 
I (l/ (l + e £ / T )). Equation (12J) can be expressed as 



dp 
8T 
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HT 



r 

'\H m + V eS , p] ~ {- , (p ~ Peg)} 



using the approximation T 

47T 3 g* T 2 
45 M„ 



(3) 

HT, where H = 



is the Hubble parameter, M p \ is the Planck 
mass and g* is the effective number of relativistic de- 
grees of freedom. Since the full Hamiltonian is com- 
plex, the equation decomposes into a coherent commu- 
tator, and an anticommutator which describes loss of 
coherence. V e s — I(V e> V^, V T , 0, 0) describes the ef- 
fects of matter on the coherent part of the neutrino 
evolution. For zero lepton asymmetry, 
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where A e = 17 and A^ T = 4.9. T = 7(r e ,r M ,r T ,0,0) 
encodes decoherence and damping due to collisions 
with the background medium, 



r Q ~ y, 



180C(3) 
7tt 4 



(5) 



with y e = 3.6, y^ r = 2.5. 

Before discussing the full numerical solution of 
Eqn. ([3]), we start with an approximate analytical solu- 
tion for guidance in understanding the physics. Briefly 
(see Appendix for more details), the problem can be 
approximately reduced to two independent equations, 
each describing the population of one of the ster- 
ile species. For each sterile neutrino, (we use v s as 
an example in the following expressions) one can ap- 
proximately use two independent oscillation channels, 
v e — » v s and — > v s . For each channel, the effective 
mixing angle in vacuum is given by 



sin 2 28 as 



4(7. 



aA > 



e,fi 



(6) 



while in-medium the mixing is suppressed according to 
the expression 



sin 2 26L 



b a (p,T) = 



sin 26 as 



(i-b a (p,T)r : 

2EV a 



A? 



(7) 
(8) 



One can then solve the evolution equation for f s , the 
phase space distribution of v s , in terms of the inter- 
play of oscillations and collisions. If f s /f a is energy- 
independent (i.e., a constant), and the mixing 9 m is 
well in the vacuum limit at the freezeout epoch, we 
find the contribution of v s to iV^P N to be 
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-2.06 x 10 3 



1.29U 



(9) 



A similar expression holds for v r , with the substitu- 
tions C/ 2 4 — > C/ 2 5 and 1714— >ms. Ultimately, the total 
contribution of the two sterile states to N c g is given by 



AN. 



BBN 
eff 



AN. 



BBN 
eff.s 



AN. 



BBN 
eflt.r 



(10) 



As expected, a sterile species is more populated at 
the time of BBN if oscillations are more efficient, i.e. 
for larger mixing (larger oscillation amplitude) and 
larger mass squared splitting relative to the active 
species, which means smaller oscillation length and 
therefore a higher probability of flavor conversion be- 
tween two successive collisions. We also stress (see 
Appendix), that AAT^ S N is larger for a larger colli- 
sion rate; indeed, collisions favor the growth of the 
sterile population toward equilibrium [38] and in the 



limit of no collisions (oscillations only), we would have 
< (sin 2 2# es + sin 2 29^ s )/2, where the right hand 
side is the sum of the average vacuum oscillation prob- 
abilities in the two channels. Note that the production 
of v s ,v r via oscillation from is more efficient, be- 
cause for the — v s system the mixing angle is less 
suppressed by the thermal potential (|V M | < |V e |). 



Let us now discuss the numerical solution. We follow 
the technique from Melchiorri et al. [llj . numerically 
evolving the neutrino densities from temperatures of 
100 MeV down to 1 MeV. To simplify the resultant set 
of differential equations, we assume a monochromatic 
neutrino energy distribution, with E u ~ 3.15T, rather 
than use the full spectrum of the neutrinos. This sim- 
plification has little effect on the density evolution [11] . 



We begin by running a loose scan across the al- 
lowed parameter space with the goal of finding ref- 
erence points that minimize N e g whilst also keeping 
m 4 and m 5 as low as possible, due to strong cosmo- 
logical constraints on the sum of the neutrino masses. 
Constraints on the masses and mixing parameters are 
from SBL data [34]. As shown in Eqn. ^}J, the con- 
tribution of each sterile neutrino to N e g is smallest 
when the two mixing matrix elements for that neu- 
trino are minimized. However, SBL constraints on 
the product of the four mixing matrix elements pre- 
vent all four elements from being small. By definition, 
m5>m4, and so it can be seen from Eqn. ^ that v r 
will have a larger contribution to A r ( ^ f BN than v s for 
comparable mixing angles. For this reason, we focus 
on minimizing the thermalization of v s . Point 1 in Ta- 
ble [T] is chosen to correspond to the minimum values of 
U e A, U^i still allowed within 2a. U e 5, f7 M 5 are chosen 
to be as small as possible while satisfying constraints 
on 4|/7 e 4 U^i U e 5 U^l- mi, are chosen to be as low 
as possible while still allowed by our choices of [Zy. 
This point does indeed lead to incomplete thermaliza- 
tion of u s , with AA r J f BN = 1.86, as shown in Figurejlj 
although this is still outside the 2a allowed range from 
Izotov and Thuan 1 1 . /V C BBN = 3.80±£j. In this re- 
gion, the degree of thermalization is quite sensitive to 
U e A, U^. As an illustration, Figure [T] also shows the 
density evolution for two additional points in param- 
eter space, labelled Points 2 and 3 in Table |TJ where 
U e A, U p_A are pushed to even lower values, outside of 
the 2a allowed region, but still within the 99 % C.L. 
allowed region. 7714, 7715 are kept fixed, and U e ^,, U^ 
are chosen to keep A\U e A U^ U e 5 U^ \ as close to the 2a 
allowed region as possible while still remaining within 
the 2a allowed region themselves. At Points 2 and 3, 
AN®fP N is safely within the Izotov and Thuan 2a al- 
lowed range. 
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shown in the Appendix, and ^ = AN^ lB . In Ta- 
ble [I] we show the effective ^2 m v f° r the three points 
we consider. 

In addition, many upper limits on the sum of the 
neutrino masses assume the standard model value of 
iV e fj, and so do not directly apply to sterile neutri- 
nos; However, a number of groups have constrained 
the N e s — ^ m v plane using various combinations of 
measurements of the CMB, Hubble constant, baryon 
acoustic oscillations, and galaxy clusters [3HH32- For 
the value of iV c ff we are interested in, detailed below, 

fn v ^ 0.7 eV is excluded at the 95% confidence level, 
which is inconsistent with the values shown in Table H] 
for all 3 points. (Tension between SBL and cosmolog- 
ical data is also discussed in Ref. [15] in the context of 
two fully thermalized neutrinos.) 

There is another, equally important factor that can 
affect the value of N e s that should be utilized when 
applying cosmological constraints: the fact that neutri- 
nos may not be fully relativistic at the time of matter- 
radiation equality. 

The standard model neutrino temperature at 
matter-radiation equality is T v = 0.55 eV [21] . [44] . 
N e g as derived from the CMB measures the relativistic 
energy density at the time of matter-radiation equality. 
A neutrino with m„ ~ O(leV) will not be entirely rel- 
ativistic at this time, and so will contribute AN e g < 1. 
Constraints on N e ff are continually getting tighter, and 
so this can be an important effect for sterile neutrinos 
towards the top of the allowed mass range. A simi- 
lar effect was considered in Ref. [IS], where the scale 
factor at matter-radiation equality was related to the 
mass and energy density of a sterile neutrino. 

The pressure density provides a convenient measure 
of how relativistic a particle is at any given tempera- 
ture, with P = p/3 for fully relativistic particles and 
P = for non-relativisitc particles. As the sterile neu- 
trinos become less relativistic, their pressure drops be- 
low p/3, and the relativistic fraction of their energy 
density drops accordingly, 




using P^/pf = 1/3 when m = 0. With this, the 
effective number of relativistic degrees of freedom can 
be expressed as 

rcl p 

= Jr- - pS^- ( 13 ) 



TABLE I. Mass, mixing parameters, results for AN e g at 
BBN and z cq , and effective mass sum for the three sample 
points discussed in the text. The derivation of the final 
three columns is discussed in Sections iHl and ITTT1 
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FIG. 1. Sterile neutrino density evolution as a fraction of 
the thermal density p cq , for the masses and mixing angles 
listed in Table [T] Dashed lines are for u a , solid lines are for 



III. PARTIAL THERMALIZATION AND 
PARTIALLY RELATIVISTIC NEUTRINOS IN 
A 3 + 2 SCENARIO: CMB CONSTRAINTS 



Both the position of the lowest peaks in CMB an- 
gular power spectrum and the damping tail at high 
multipole moments, I, are sensitive to the redshift at 
matter-radiation equality, and hence to the energy den- 
sity of relativistic neutrinos at that time. The latter 
is usually expressed as a constraint on the sum of the 
neutrino masses, 



^m v = 94eV(^ m /i 2 ), 



(11) 



where VI is the density as a fraction of the critical den- 
sity of the Universe, fl = p/p c , and f2„.m is the neu- 
trino contribution to the matter abundance f2 m . It 
is important to note that in Eqn. ( |11[ ) it is assumed 
that each species is fully thermalized, by assuming 



that for each species, 



-rel 



77i„rt* h , with 



,th 



from 



Eqn. (14). Constraints on the sterile neutrino mass are 
really constraints on the product m v n v , and if a ster- 
ile neutrino does not undergo full thermalization, then 
it contributes m e ff = m v to constraints on ^ m v . 
For the partially thermalized m„ — 0.6 eV neutrino 
we consider in Table [T] the phase space distribution 
is approximately a scaled Fermi-Dirac distribution as 



The number, pressure and energy densities are given 
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by the standard formulae, 



Pu = 



9 
2tt 2 

2^2 / ''I'TE-t'-U'^- 



9 
2tt 2 



dpp 2 fu(p), 
P A 



d P Ep 2 U(p), 



(14) 
(15) 
(16) 



where g counts the number of helicity states, f v (p) 
is the neutrino phase-space distribution, and p = \p\. 
When the neutrinos are in thermal equilibrium, they 



follow a Fermi-Dirac distribution, 

Up) ' 



1 + exp[f ] ' 



(17) 



when the chemical potential is zero. After freezeout at 
Tp ~ 2 MeV, the co-moving number density must be 
conserved: 

n v {p,T) = ) n„(p F ,T F ) 

V a J 2-k z J 1 + exp[^] 

where a is the scale factor, and subscript F denotes 
the value at freezeout. Neutrino momentum redshifts 
as p — (^f 1 ) PFi and so the neutrino number density 
after freezeout is 




Comparing this with Eqn. (14) and using T = T /a, 
we have 



U(p) 



1 + cxp 



V 



T, 



(20) 



after freezeout. This reduces to the standard expres- 
sion for both relativistic and non-relativistic particles. 



iVcff can then be found using Eqn. ( 13 ) with Eqns. ( 15 
20 ) . Since the co- moving number density is conserved 



from the time when the neutrinos were entirely rela- 
tivistic, it follows that the number density dn/dp must 
be independent of mass. Thus the total neutrino en- 
ergy density dp/ dp — \/p 2 + m 2 dn/dp will be larger 
for particles with larger mass. This is compensated 
for by the fact that the pressure is smaller for massive 
neutrinos, leading massive sterile neutrinos to have a 
smaller contribution to N e s as shown in Figs. [2] and [3j 
We use this relation to determine the relevant val- 
ues of iV e ff at matter-radiation equilibrium for the 
3 + 2 sterile neutrino models discussed earlier to ex- 
plore whether they may be consistent with both SBL 
and cosmological bounds, and report the results in Ta- 
ble 1 @B]. The use of Eqn. (13) to determine N e g 



requires knowledge of the neutrino phase-space distri- 



bution at decoupling, and we demonstrate in the ap- 
pendix that our three points satisfy the conditions re- 
quired for the phase-space distribution to be approxi- 
mated by a Fermi-Dirac distribution scaled by a con- 
stant. With this approximation, and using Eqn. (13) 
with Eqns. ( 15|20[ ), the contribution to iV c 
AiV eff = (1.68, 1.47, 1.25) at Points 1, 2, and 3 respec- 
tively. This leads to some easing of the tension between 
SBL data and CMB constraints on N B g. 



IV. DISCUSSION AND CONCLUSIONS 

We have discussed the cosmological generation and 
evolution of a population of two sterile neutrinos, with 
masses and mixings motivated by the recent SBL data. 
Specifically, we calculated the contribution of these ex- 
tra neutrino species to N e g at BBN and CMB epochs. 
We focused on the region of the parameter space where 
the sterile neutrinos are produced with less than ther- 
mal abundance (AN c g < 2), so that the tension with 
BBN and CMB measurements is eased, compared with 
the case of two fully thermalized species. We find 
points at the limit of the region of parameter space al- 
lowed by the SBL data where the heaviest sterile state 
is fully thermalized, while the second is produced with 
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FIG. 2. Contribution of one thermalized massive sterile 
neutrino to N e g at the time of matter-radiation equality. 
If the sterile neutrino has an approximately Fermi-Dirac 
distribution, this is equivalent to AN^ / AJV < ^ j P N , i.e. if the 
sterile neutrino is not fully thermalized and ANeg < 1, 
then AN^ will be reduced accordingly. 




FIG. 3. Contribution of one massive sterile neutrino to iV e ff 
as a function of the equivalent temperature of a massless 
neutrino, T„. At matter-radiation equality, T v = 0.55 eV 



abundance as low as ~ 40% of the thermal abundance. 

Whilst it is possible - with the maximum suppression 
of N e g due to partial thermalization - to find points 
in parameter space marginally compatible with BBN 
constraints (N^ N < 4.6 at 2a), the tension with BBN 
data overall remains. 

Interestingly, if SBL-favored sterile neutrinos really 
are the origin of Neg > 3, we expect their contribu- 
tion to N e g at z eq - relevant for CMB constraints - 
to be lower than that at BBN epoch, due to their be- 
ing only moderately (partially) relativistic at z eq . In 
principle, this feature would allow us to distinguish 
the sterile neutrino hypothesis from other possible ori- 
gins of an excess of radiation, however the difference 
AtBBN _ AfCMB ig on the order of 1Q % or less. It is not 



clear whether foreseeable probes of N e g at different 
epochs will ever reach this level of precision. We also 
note that while the mass-induced suppression works 
to ease the tension with the CMB data somewhat, it 
comes with a price: the ~ 1 eV masses of the sterile 
states would increase the sum of the neutrino masses 
to m v > 1 eV, which is disfavored by CMB bounds 
on this quantity. 

Summing up, we find that even with the suppression 
effects due to partial thermalization and partially rel- 
ativistic masses, two additional sterile neutrinos in a 
mass range that might explain SBL neutrino data ap- 
pear to be inconsistent with cosmological bounds com- 
ing from BBN and CMB measurements. 

It has recently been questioned whether the SBL 
data from MiniBooNE actually favor 2 sterile neutrinos 
|31j . If this requirement is relaxed then the results we 
derive will be particularly relevant to constrain mod- 
els with one extra neutrino. Alternatively, new physics 
that might resolve these inconsistencies include alter- 
ing microphysics or altering cosmology. An example 
of the former includes introducing a lepton asymmetry 
in addition to the two sterile neutrinos [4"TH4"5] . For 
example, it has recently been shown that large lepton 
asymmetries do lead to incomplete thermalization in 
a scenario with one sterile neutrino [7]. As an exam- 
ple of the latter, Hamann et al. [115] find that mod- 
els with one fully thermalized eV-scale sterile neutrino 
and additional massless degrees of freedom can pro- 
vide a better fit to a wide range of cosmological data 
than standard ACDM, if the dark energy equation of 
state parameter is free to be w < — 1. Nevertheless 
this requires a contribution to N c g from massless ster- 
ile states of AN e g > 1, and is thus in tension with BBN 
constraints on Y p , unless the standard model neutrinos 
have a non-zero chemical potential. 

The need to consider such exotic possibilities to pos- 
sibly obviate the bounds we derive here demonstrates, 
once again, the important utility of cosmological obser- 
vations on constraining fundamental neutrino particle 
physics. 



ACKNOWLEDGMENTS 

LMK and TDJ acknowledges support from the DOE 
for this work; CL and TDJ acknowledge the support 
of the NSF under Grant No. PHY-0854827. 



Appendix A: Analytical Approximation 



AN*™ 



Let us derive th e a pproximate analytical result for 
, Eqs. @-jl0f. 
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Oscillation amplitude in vacuum and in 
medium 



We consider a basis of 5 neutrino flavor states (v a , 
a = e,[i,T,s,r, where v s ,v r are the sterile states. 
These are related to the mass eigenstates by the mix- 
ing matrix U, v a = £) i=1 5 U m v l , where U s ± ~ U r5 ~ 1 
and U is taken to be real for simplicity. As stated in the 
main text, we use a number of assumptions to simplify 
the problem: (i) U s5 = U ri = (ii) f/ T 4 = U t5 = 0, 
(iii) mass hierarchy > m 4 3> rrij. 

Let us find the amplitude of active-sterile oscillations 
in vacuum, by calculating P[y a —¥v s ) (for defmiteness; 
analogous results hold for v r ), with a = e, fi. The chan- 
nel v T —> v s has zero probability due to assumption (ii) . 
We use the standard notation Am 2 , = m 2 — m 2 . 

In all generality, the standard oscillation formalism 
gives 

/ Am 2 - \ 

P{v a -¥u B ) = -4 V U ni U,,iU nA U S 4 sin 2 [ -^-t 



/ ^ ^ ai ^ si ^ aj U sj 1 
i,j,i>3 



4E 



(Al) 

Using assumption (i) and neglecting the lowest oscilla- 
tion frequencies (assumption (iii)), we get: 

P{v a -*v a )~ -4U a4 U s4 ^2 U aiU sl sin 2 (u s t) , 



i=l,2,3 

4[/ 2 4 t/ 2 4 sin 2 ( w>t j 



(A2) 



where u> s = rn^/(4E), and the last expression is ob- 
tained using the unitarity of the mixing matrix. Eq. 
( A2 ) has the same form as the classic 2-neutrino oscil- 



lation probability, with effective mixing 



sin 2 20„, = AUlAJ, 



a4 u s4 



4U: 



a 4 



(A3) 



An analogous result is obtained for v r , with u r = 
m^/(4E) and sin 2 29 ar ~ 4C/ 2 5 . It is immediate to 
verify that, under the same assumptions as above, 
P(y a —±v r ) — P(i/ r —¥i/ a ) and P(i/ a —^v s ) = P(v s ^v a ). 

Due to the thermal refraction potential, the effective, 
two-neutrino oscillation amplitude is suppressed - for 
both neutrinos and antineutrinos - as follows (see, e.g., 



sin 2 29 m 
b a {p,T) 



sin 2 29 a 



(i-b a (p,T)r 

2EV a 



Am 2 



(A4) 
(A5) 



where V a is as in Eqn. Q and we used cos26* as ~ 1, 
for convenience in the calculations that follow. The 



expression ( A5 1 is valid for a CP-symmetric neutrino 



gas; the more complicated case with a lepton asymme- 
try will not be discussed here. 



2. Flavor evolution equation and its solution 

Let us now consider the production of the sterile 
neutrino v Sl using an effective two-neutrino system, 
v a — v s with the oscillation frequency and amplitude 
as outlined above. Let f s and f a be the phase space 
distributions of v s and of one of the active species, and 
let p be the neutrino momentum. We start with the 
evolution equation (see e.g., Foot and Volkas [5U] and 
Dodelson and Widrow EH), 



li- HE -m)^ EM 



sin 2 29 m (E,t) 



r a {E,t) 



(f a (E,t) - f s (E,t)) (A6) 



with r Q being the collision rate, Eqn. Eq. (A6) 



is valid when the neutrino oscillation length is much 
shorter that the neutrino mean free path, so that 
the effect of oscillations between two collisions is de- 
scribed by the averaged in-medium oscillation proba- 
bility (P) = sin 2 26> m /2. We have verified that this 
is always the case for our parameters of interest (52) . 
We take f a = f a {p/T) to be a Fermi-Dirac distribu- 
tion, but the derivation in this section holds for any 
function of p/T. 

Eqn. (A6) can be simplified using [ST] 



d d 
dt~ HE dE 



= -HT 



d_ 
df 



y = p/T(A7) 



and defining 



o5/4 <p3 

— — V^U(3)G F — 
irMw m.4 



3.53 x 10 



so that —b a (p, T) 
equation 




(A8) 



T y /eV" 
MeV J Vm 4/ 

Thus, we have the new 



HT 



df s _ sin 2 29 as T a (y,x) 
dT 



2(1 



x 2 y 2 ) 2 



(f a (y)-fs(x,y)), (A9) 



where y and x should be treated as independent vari- 
ables, and the differential equation should be solved 
with respect to x, with y fixed. 

Changing the differentiation variable from T to x 
[5"T] , and neglecting the temperature dependence of g* , 
one finds a solution of the form 



1 



k 

fa 

to/: 



exp 



TO4 

13V5 5 
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G 3 J 2 M pl M w 
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1.05 x 10~ d eV , (A10) 



where g* = 10 has been used. As discussed in Ref. [ST] , 
if the lower integration limit is small, i.e. xy <C 1 



8 



(—b a (p,T) -C 1) at freezeout, we can replace it with 0, 
thus solving the integral exactly, with 7r/4 as the result. 
In this limit, f s /f a is independent of y [51], meaning 
that f s has the same spectral shape as f a , and the 
two only differ by an overall factor. For the freezeout 
temperature T„ ~ 2 MeV, and for p = 3.15T, we get 
xy ~ 4 x 10~ 3 m 4 /eV, so this condition holds for the 
range of masses of interest here. 

The final result for the ratio f s / f a is then 



k 

fa 



exp 



7T m 4 2 

4 TTlf 



JJa 



(All) 



immediate: 



fa 



1 — exp 




exp 
exp 



7r m 4 

1 /// ; \ " e " N '71 

-2.06 x 10 3 



U 



1.29L/2 4 ) 



(A12) 



A similar formula holds for the abundance of i/ r , 
upon replacement of index: 4—^5. Note that the 
result in Eq. (A12| can be rewritten in terms of a 
single, effective. 



sin 2 26 eS = 4(Ui 



— v s system, with mixing angle 
1.29C/ M 2 4 ). 



In the limit f s /fa *C 1, an expansion of the exponential 
recovers the result in [TTJ [ST] . 

For our case, where two oscillation channels are 
present, v s «-> v e and v s •<->• v^, the generalization is 



We find that our analytic solution, Eqn. (A12), gives 



results around 10% lower than our numeric solution at 
Points 1,2,3 in Table|T] The main source of this discrep- 
ancy is g* , which is kept fixed in the analytic solution, 
while its full temperature dependence is included in 
our numeric results. When g* is kept fixed in both 
calculations, results match to within 5%. 
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